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Abstract. Consider the initial-boundary value problem for the 2-speed 
Carleman model of the Boltzmann equation of the kinetic theory of gases 
(see [T. Carleman, "Problemes mathematiques dans la theorie cinetique 
des gaz", Almqvist-Wiksells, Uppsala, 1957]), set in some bounded inter- 
val with boundary conditions prescribing the density of particles enter- 
ing the interval. Under the usual parabolic scaling, a nonlinear diffusion 
limit is established for this problem. In fact, the techniques presented 
here allow treating generalizations of the Carleman system where the 
collision frequency is proportional to the a-th power of the macroscopic 
density, with a G [— 1, 1]. 



In the 1930's, Carleman proposed a model [1] describing the time evolu- 
tion of a mono dimensional gas composed of two species of particles that move 
at a constant speed c > in the j;-direction. The number density at time 
t and position x of particles moving at speed +c is denoted by u = u{x, t) 
while that of particles moving at speed — c is denoted v = v{x, t). Carleman's 
system is 



A very natural extension of the model, called the generalized Carleman 
model, involves a collision frequency that is proportional to some power of 
the macroscopic density p = u + v, as follows: 



For a = 1 the original Carleman system is recovered, whereas a = gives 
another remarkable system, known as the Goldstein- Taylor model [Sj I21j . 
that can be reduced to a damped wave equation (the telegrapher's equation). 
In the latter case, there is an explicit representation of the solution in terms 
of the standard Poisson: see \12j, p- 56 for this formula, originally found by 



1. Carleman models and their diffusion limits 




dtu + cdxU = {u + v){v — u) , 
dtv — cdxV = {u + v){u — v) . 




dtu + cdxU = (n + v)°'{v — u) , 
dtv — cdxV = {u + v)'^{u — v) . 



M. Kac. 
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Other instances of kinetic models with a colhsion frequency whose de- 
pendence on the macroscopic density is other than hnear can be found in 
Radiative Transfer: see for instance |13j . 

The only nontrivial hydrodynamic limits of the Carleman model, general- 
ized or not, are diffusion limits — linear or nonlinear. Indeed, local equilibria 
for those models are those densities for which (n — v){u + u)" = 0, implying 
u = V. Hence all local equilibria for those models have mean velocity 0: in 
the language of the kinetic theory of gases, analogues for the system ([T]) of 
local Maxwellians have mean velocity 0. Hence the limiting dynamics of the 
system can be observed only on a longer time scale that corresponds with a 
diffusion equation. 

After setting c = 1 without loss of generality, we consider the limit as 
e ^ 0+ of the following scaled version of ([2]): 

e'^dtUe + ed^Ue = (n, + - u,) , 

e'^dtv^ - edxVe = (n^ + Ve)"{ue - v^) . 
Defining the macroscopic mass density pe and the current by 

(4) Pe = Ue+V^, J, = , 



we put the system ([3]) in the form 

dtPe + dxje = , 

e^dtje + d^Pe = -2p"je ■ 

Then, we prove that the term e^dtje is negligeable in the vanishing e limit, 
and show that the limiting density is governed by the following nonlinear 
diffusion equation 

(6) dtp = \dxx ( 

for a E [—1, 1), while the case a = 1 leads to 

(7) dtp = \dxx^'n. p . 

Several result on this problem have been obtained over the last thirty 
years. 

In particular, existence and uniqueness for the solution of the initial- 
boundary value problem of ([1]) have been proven by Fitzgibbon [3] — see 
also [7] for the same problem on the infinite line, and [16], [T7] and EO] for 
more information concerning the large time behavior of the solutionqj. 

Several authors (among others Kurtz [8], McKean [12], Fitzgibbon [3j, 
Pulvirenti and Toscani [15], P.-L. Lions and Toscani [9], Marcati and Rubino 
[TT] , Donatelli and Marcati [2] , Salvarani and Vazquez [19] ) have considered 
the relaxation problem for the Carleman system or generalizations thereof. 



^We are grateful to L. Tartar for indicating the references jTj and [16], [17] . 
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All these papers on the diffusion limit for kinetic models deal with the 
initial value problem on the infinite line, or with the initial-boundary value 
problem with specular or periodic conditions at the boundary. 

Establishing the diffusion limit for Carleman type models (generalized or 
not) in a bounded domain with non homogeneous boundary conditions is 
more difficult. Indeed, mimicking the proofs of the diffusion limit in the 
cases mentioned above fails to produce a uniform bound on the current je 
in the case of a nonhomogeneous boundary value problem. 

To the best of our knowledge, the only available convergence proof for an 
initial-boundary value problem with very general boundary conditions at the 
time of this writing is in [T8] for the Goldstein- Taylor model. Unfortunately, 
the proof in [18] uses extensively the linear nature of the problem coming 
from the assumption a = 0, and does not seem to be extendible to nonlinear 
models. 

The present work establishes the diffusion limit for all Carleman systems 
([2|) with Q £ [—1, 1], in a bounded domain 0, = (0, 1), with boundary condi- 
tions imposing the density of particles entering the domain fl 

u,{t,0) =ip-{t), t>0, 

v,{t,l) =ip+{t), t>0. 

We have restricted our analysis to the case \a\ < 1, as it leads to a unified 
treatment based on the dissipative nature of the problem. 

The new ingredient in the present paper is a uniform estimate bearing on 
some notion of relative entropy of the solution (u^, fe) with respect to a well- 
chosen profile that satisfies the boundary conditions expected to hold in the 
vanishing e limit. By failing to introducing this profile, previous attempts 
to establishing the nonlinear diffusion limit in the case of nonhomogeneous 
data fell short of obtaining uniform bounds on the current j^, except in the 
particular case f~^{t) = ^~{t) = Const. 

This idea of using the relative entropy with respect to some adequate pro- 
file in order to match nontrivial boundary conditions can be used on more 
complicated models. In the introduction above, we already mentioned that 
the generalized Carleman equations are somewhat analogous to more com- 
plicated models appearing in Radiative Transfer, for which the nonlinear 
diffusion approximation is known under the name of "Rosseland approxima- 
tion" . We shall apply the method presented here to these more complicated 
models in a forthcoming paper 

2. Main results 

Before stating the convergence theorem that is the main result in this 
paper, we recall some background on the Carleman systems. 

2.1. Existence and uniqueness theory for Carleman systems. The 

class of initial and boundary data for the Carleman systems considered in 
this paper, which we henceforth call "admissible data", is defined below. 
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Definition 2.1. Consider the system @ with \a\ < 1 posed for {t,x) G 
(0,T) X (0,1), with boundary conditions (0) and initial condition 

ne(0,x) = , 0<x<l, 

VeiO,x) =v'''{x) , 0<x<l. 

The initial data and the boundary data {ip~,(p^) are said to be 

admissible on the time interval (0, T) if and only if 

(1) ip^ £ W^'°°{[0,T]) and ip^ > on {0,T), while 

(2) u*" and v'"" £ (0,1) and u^", > a.e. in (0,1). 

The existence and uniqueness of a nonnegative solution for the original 
Carleman model is well-known. The following result is a standard general- 
ization of Fitzgibbon's in [4j: 

Theorem 2.2. Consider the generalized Carleman model ([H) for \a\ < 1 in 

the domain Q = (0, 1) with boundary condition and initial condition 
Assume that the initial and boundary data (u^^,v^"') and ip^ are admissible. 
Then there exists a unique nonnegative generalized solution (n, v) of ^ in 
C([0,r];Li(0,l) xLi(0,l)). 

Proof. It is well known that the unbounded operator 

Ba{u, v) = {—dxU + {u + v)'^{v — u),dxV + (n + u)"(n — v)) 

on ^-"^(0, 1) X L^{0, 1) with domain 

V{Ba) = {(n, v) e ^^'^(0, 1) X W^'^iO, 1) I u{0) = and v(l) = 0} 

is dissipative (see, for example, [HI US]). The existence and uniqueness of 
the solution {u,v) easily follows from the same method as in [1]. □ 

2.2. Uniqueness theory for the nonhnear diffusion equation. Next 
we consider the target, nonlinear diffusion equation ([6]) (or ([7]) when a = 1). 

Definition 2.3. Consider the nonlinear diffusion equation (0) or Q when 
a = 1 with Dirichlet boundary conditions 

(10) p(t,0) = 2(^-, p(t,l) = 2(^+, tG(o,r), 
and initial condition 

(11) p{0,x)= po{x), XG(0,1). 

We call these data admissible if and only if 

(1) for a G [—1,1), po > belongs to L^((0, 1)) while ip^ > is in 

Tyi'~([o,r]); 

(2) for a = 1, same assumptions except that ip^ > 0. 

The notion of admissible weak solution of the nonlinear diffusion equations 
([6]) or ([7]) being less obvious than in the case of the Carleman model, we 
recall it below. 
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Definition 2.4. An admissible weak solution of ^ or Q with Dirichlet 
boundary conditions 

(12) p{t,0)=2ip-, pit,l)=2ip+, tG(0,r) 
and initial condition 

(13) p{0,x)=po{x), xG(0,l) 

for admissible data {po,ip^) is a nonnegative element of L'^{{0,T) x (0,1)) 
such that dc.p^-" £ L^{{0, T) x (0, 1)) if a / 1 while Inp G ^L((0> ^) x (0> 1)) 
with dx\n.p G L'^{{<d,T) x (0,1)) if a = 1, that satisfies for a.e. t G 
(0,T), together with the identities 

Z'^*'?^ ~ 2{i-a) ^^P^~°'^^^) ^)^2;cit + / poix)(l){x, 0)dxdt = , 



JO 



for all (j) G C([0,r] x [0, 1]) n H^{[0,T] x [0, 1]) vanishing in x = 0, x = I 
and t = T if a ^ 1, whereas, for a = 1, 



(15) / / [pdt<j) — ■^dxlnpdx(l)) {t,x)dxdt + / po{x)(l){x,0)dxdt = . 
Jo Jo Jo 

Notice that, in what fohows, we only need to know that the boundary 
value problem for ([6|) or ([7|) has at most one admissible solution, as the 
existence of such solutions will result from that of a solution of the general- 
ized Carleman model in the hydrodynamic limit. The proof of the following 
lemma is classical: 

Lemma 2.5. For each set of admissible data (pQ,ip^), the nonlinear diffu- 
sion equation @ for — l<a<l orQ if a = 1 has at most one admissible 
solution. 

2.3. The convergence result. The main result in the present paper is the 
following convergence theorem. 

Theorem 2.6. Let (ti™, c/?^) be admissible data for the generalized Car- 
leman system on the time interval [0, T]. For each e > 0, let (ne,Ue) be the 
solution of the scaled Carleman system with initial condition Q and 
boundary data ([2|. 

Then, in the limit as e — > 0, the macroscopic density 

p, = u, + v,^p ^nLL([0,^] X [0,1]) 

where p is the generalized solution of ^ if a £ [—1, 1) or o/ Q if a = 1, 
with initial and boundary conditions 

p(x,0) ='U*"(x) +i'*"(x), xG(0,l), 
p(t,0)=2vp-, p(t,l)=299+, iG(0,r). 
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3. Uniform bounds on the density and current 

As mentioned in the introduction, the uniform (in e) bounds on the den- 
sity and current are consequences of the equation that is satisfied by some 
notion of relative entropy, which we shall introduce below. 

From now on, we set 

(16) V = max(||(^+||^yi,oo(o,T), \W~\\w^^^{o,T)) 
and 

(17) ip.„, = mi{ip^{f) I < t < T} > . 

We begin with some elementary background on convex functions. 

3.1. A family of convex functions. Let ^ : R+ — > R be a convex func- 
tion that is on R^, and that satisfies 

M < 

and 

> +00 as y ^ +00 . 

y 

Define the Legendre dual (j)* of (j) by 

cp*iO = ^M^y-Hy))- 

y>0 

This definition clearly implies that 

< Hy) + <f>*iO for each y,^eR+. 
This inequality is a generalization of Young's classical inequality for 

cj>{y) = ^ and <P*iO = ^ 
p p' 

where p and p' are dual exponents in the sense of Holder's inequality: 

p G [1, oo) and p' = . 

p—l 

More specifically, we shall use the family of convex functions below: 

a) for P G [—1, 1), the function (?!)^ is given by 

= 2^y'~^ , y>0; 

b) for /3 = 1, the function (^i is given by 

0i(y) = ylogy and 0i(O) = 0, y>0. 
Straightforward computations show that 

-, a 2-/3 

while 

for each ^ > and P G [-1, 1). 
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With V as in (fT6]) . for each /? G [—1,1], we set 

(18) ci>^(y) = 0;3(y) + '/'M'/';3H + i), y>o. 

Young's inequahty entails 

$/3(2/) > (0;3(z^) + l)y, y>0, 

and since = <^^^ one has 

(19) $^(y) - ^'p{v)y >y, /? € [-1, 1] , y > . 

3.2. The relative entropy and entropy production rate. Our bounds 
on the Carleman system involve the notion of relative entropy with respect 
to some suitable profile that satisfies the same boundary conditions as the 
solution of the limiting diffusion equation. There are many possible ways of 
choosing this profile. When the boundary data are constant, a logical choice 
would be to pick the stationary solution of the diffusion equation with those 
boundary data. In the case of time dependent boundary conditions, this 
choice is less natural, and we simply define this profile to be the convex 
combination of boundary data: 

f{t,x) := {l-x)ip-{t)+xip+{t). 

Given $ : R_|_ ^ R, a C'^ convex function, we define the relative entropy 
of the 2-velocity density {u,v) = {u{x),v{x)) E (R+)^ to be 

(20) Hp[u,v\f]= [\'^f,{u) + '^(,{v)-2<^f,{f)-^'p{f){u + v-2f)]{x)dx, 

Jo 

and the entropy production rate as 

(21) Pp[u,v]= f {<!>'p{u) -<^'fj{v)){u + v)°'{u-v){x)dx>0. 

Jo 

Assuming that (ue,Ve) is a solution of the scaled Carleman problem 
we next determine the equation governing the evolution of Hf^[ue,Ve\f]. Mul- 
tiplying each side of the first equation in ^ by $^(tie) and each side of the 
second equation in ^ by ^'^(v^), one finds, upon adding both resulting 
equalities and integrating in x over [0, 1]: 

d 1 fi 

— / {<^>p{u,) + <^>p{v,)){t,x)dx + - d^{<^p{u,) -^i3iv,)){t,x)dx 
"f^ Jo e Jq 

1 „ r 
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It is actually more convenient to rearrange this equality as follows: 
d 

- {<S>p{u,) + ^fs{v,)-<^>'^if){u, + v,))it,x)dx 

1 

+- 

(22) 



+- / {<l>p{u,) - <i>p{v,) - ^'p{f){u, - V,)) {t,x)dx 



1 1 

-^Pp[Ue,Ve\ / dx^l3if){Ue - Ve){t,x)dx 

+ / dt^i3{f){u, + v,){t,x)dx. 



Observe that 
1 

- / dx {^^{Ue) - '^p{Ve) - <^p{f){u, - V,)) {t,x)dx 

f Jo 

= l[^f,{u,)-^^{f)-<^'^{f){u,-f)]l 

(23) --^[^^{v,)-^fs{f)-^'^{f){v,-f)]l 
= i (cl>^(n,) - <^p{f) - ^'fs{f){u, - /)) {t, 1) 

+1 _ ci,^(/) _ <!>'^{f){v, - /)) [t, 0) 

since 

n,(t,0) = /(i,0) and v,{t,l) = fit,l) . 

By convexity of both terms in the last right hand side of (|23p are non- 
negative; hence 

1 

(24) - / dx i^fsiu,) - <^^{v,) - <^'p{f){u, - V,)) {t, x)dx > . 
e Jo 

Next, we formulate the equality (j22p in terms of the relative entropy 
H^lu^^Velf] as follows: 

4:Hl3[Ue,Ve\f] + \Pf3[Ue,Ve] < / d^^^ i3{f){Ue - Ve){t, x)dx 

(25) ^ ' ° 1 

+ / dt^'p{f){u, + v,){t,x)dx+ j ^'l,{f)dtf\t,x)dx. 



We estimate the first term on the right hand side of the inequality (j25p in 
terms of the current je = -{u^ — v^), as follows: for each 7 G (0, 1), one has 
(26) 
1 



d,'^pif){u,-v,){t,x)dx<^ Ut,xfdx + ^\\d.,(^p{f){t,x) 

Then, the inequality (fT9]) satisfied by <^/3 implies that 

(27) ^p{y) - ^'p{f{t, x))y > ^p{y) - ^'p{v)y > y 
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for each y > 0. Indeed 

f{t,x) < max(||v9"||ioo(o,T), \\^'^\\l°°{o,t)) 

so that 

since is nondecreasing, being convex. Because of (|27p . one has 
(28) 

dt<^p{f){Ue + Ve){t,x)dx 

Jo 

< \\dt^'p{f)\\L^{Hp[u,,v,\f] + 2\\<^>pif) - 1>^(/)/|1l-) . 







With this estimate for the second term on the right hand side of (j25p . we 

recast this inequahty as 

(29) 

j^Hf,[u,,v,\f] + ^PpK,v,] <l j,{t,xfdx+\\dt^'p{f)\\L^Hp[u,,v,\f] 

+^ii9x$/3(/)iii- + ii^;^(m/2|iL- 

+2\\dt^'p{f)\\L^\\^p{f)-^'p{f)f\\L^. 

3.3. Current estimate. Our first objective is an L^^ bound on the current 

je = \{Ue - Ve). 

Start from the relative entropy inequahty (j29p with the particular choice 
(3 = a (we recall that a is the exponent in the nonlinearity of the Carleman 
system 

We are going to show that, for 7 small enough, the term 

^ / je(t, x)^d2; is dominated by -2 -Pq [tie, fe] • 
JO £ 



Indeed, if a = 1 



1 1 /^^ 

-lTPa[Ue,V^] = ^ {ul- vl){\ogUe - logVe)dx 

^ -2, , slog Me - logUe 
Ue - Ve 

> / Je > / jfdx. 

Jo max(Ue,Ue) Jo 



in view of the elementary inequality 



— > 7 — TT , for each a,b > . 

a — max(a, 0) 
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For a G [0, 1) 



1 1 /' 

1 H^~°' — 

f^{u,+v,r^ ^—dx 

Jo max(ue,u,)" Jo 

Finally, for a £ [—1,0), we separate (0,1) into the region where Ue(t,x) > 
v^{t,x) and its complement, where Vf:{t,x) > u^{t,x). 



1 1 /"^ 

-^Pa[Ue,Ve] = ^ {Ue- Ve){Ue + Ve)" {ul~" - vl~'^)dx 

Jo 



f^{u, + v,r^ — 



1-a _ 



The first integral is estimated as follows: 



1— a _ a /" _ „,l-a 

+ ^J" ^- dx > / jl{2uX ^ dx 

„ l+|ci| l+|a| 

= 2 / J, -,,11 r-^dx 



> 2" / j.^dx 

The second integral is estimated similarly. In all cases, 



(30) ^Pa[Ue,Ve]>Ca [ 



Jedx 



with 

Ci = 1 while Ca = 1 - a if < a < 1 and C„ = 2" if - 1 < a < . 

Proposition 3.1. Let {ue,Ve) be a solution of the scaled, generalized Car- 
leman system (0) with admissible initial and boundary data. Then there 
exists a positive constant J = J{u,(pm,Oi,T,u^^ ,v'^^) > such that the cur- 
rent je = -{ue — Ve) Satisfies 



rT rl 

■2 



{t, x)dxdt < J 



Jo 



for each e > 0. 



DIFFUSION LIMIT FOR GENERALIZED CARLEMAN 11 

Proof. Start from (f29|) . and pick 7 = Ca/2. Upon integrating both sides of 
(j29p over the time interval [0,t], one obtains 

H^[u„v,\f]{t) + ^ [ [ Us,xfdxds 
Jo Jo 

< \\dMf)\\L^ [ H^K,vMs)ds 

Jo 

+ Mo{u,a,^n,)t + Ha[u''',v'^\f], 

where 

+2\\dMmL^\\^M)-Kif)f\\L^- 

Gronwah's inequahty imphes that 







from which we deduce that 



^0 



3C*q; 

□ 



3.4. Density estimate. With the current estimate at our disposal, we next 
obtain an L^^^(R+; L^) bound on the macroscopic density = + v^. 

To do this, we apply ([29|) with /3 = 0; here, the entropy production rate 
is useless so that we actually discard it from the left-hand side: 



-i^O [n„ t^e I /] < i J. (t , + 1 1 a* CD^ (/) 1 1 ioo i^o [ne , I /] + Ml (i/, 7) 



d_ 
di' 
where 

Mi(z.,7) > ^wd^fwl^ + Wdtfh^ + (1 + i^fWdtfh^ 

>^||9,$o(/)||i^ + ||^o(/Wlli- 

+ 2\\dt<^>'oif)\\L^\mf) - '^oif)f\\L^ ■ 

Integrating the relative entropy inequality above over the time interval [0, t], 
we obtain 

Ho[u,,v,\f]it) < i7oK",t;'"|/] + piu,^m,a,t,u''',v'^) + M,iu,^)t 

+ \\dt^'o{f)\\L^ [ Ho[u,,vMs)ds. 



Gronwall's inequality implies that 

< {Ho[u'^,v'^\f] + lJii^,ipm,a,t,u'^,vn + Miii^,j)t) e*ll^**o(/)ll 
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On the other hand, 

Ho[u„v,\m) = i - f{t)\\l, + \\v,{t) - /(t)lli,) . 

Summarizing the discussion above, we have proved the following density 
bound. 

Proposition 3.2. Let {u^,Ve) be a solution of the scaled, generalized Car- 
leman system Q with admissible initial and boundary data. Then there 
exists a positive constant K = K{u^ipm,Oi,T,u^^ ^v^"^) > such that the 
macroscopic density = Ue + satisfies 



[ [ pl{t,x)dxdt < K 
Jo Jo 



for each e > 0. 



3.5. An entropy production estimate. We shall conclude this section 
with a further nonlinear estimate that eventually controls some flux in the 
diffusion limit. 

We distinguish the cases a € [0, 1] and a G [—1,0). 

If a G [0, 1], we use 1^ with /? = ±1. 

Indeed, the entropy production rate in the case /? = — 1 is 

P-i[u„v,] = f {u, + v,r{u, - - <^'_i{v,))dx 

K + v,nu, - v,){u', - vl)dx = C p\+^{ej,fdx . 



Hence the relative entropy inequality (I29p with f3 = —1 becomes 
d 



dt 



H^i[u,,v,\f]+ f p\-^''jUx<l f 
Jo Jo 



Ut,xydx + M2{u,-f) 



+ \\dt^'-i{f)\\L^H.i[u,,v,\f] 

with 

M2(i^,7) > 2||9tci>'_^(/)||^^||ci>_,(/) _ ci>'_^(/)/||^^ 

+ ^\\d,<^_,if)\\l^ + \\.^'L,{mf'\\L^. 

Integrating the relative entropy inequality above in the time variable, and 
applying Proposition 13. 11 we obtain 

H-i[u,,v,\f]it)+ f C p\-'''fjx<\\dt^'_^{f)\\L^ [' H-i[u,,vMs)ds 
Jo Jo Jo 

so that, by Gronwall's inequality, we arrive at an estimate of the form 



JO 



< 
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Doing the same with (3 = 1 leads to 

Pi[u^,v^] = / {ue +Ve)°'{ue - u,)(ln(ti,) - ln{v^))dx 
Jo 



Jo Ue - Ve Jo 

The estimate above rests on the mean value theorem: 

a-6)(lna-ln6) > — —^>- , a, 6 > . 

max(a, b) a + b 

Integrating in t the relative entropy inequality (|29p with (3 = 1 and applying 
the Gronwall inequality as above gives 




JO 



where 

M3{U, ^n., 7) > 2\\dMf)\\L^ W^iif) - <I>'i(/)/||l» 

+ ^\\dMf)\\i^ + \mmf\\L^. 

If on the other hand a S [—1,0), we use (p9|) with (3 = \a\. Notice that 
$1^1 is decreasing on H*^: indeed 

$f„|(z) = (l-|a|)z-H ifae(_i,o), 
$f„,(z) = z-^ ifa = -l. 

Hence, for each a,b > 0, the mean value theorem implies that, for some 

(0,1) 

(a - 6) (cDj^i (a) - $|,| (5)) = (a - b)Hl^{{l -e)a + eb)>{a- b)Hl^ (a + b). 
Therefore 

P\a\[Ue,Ve] = (u, + f ,)" ^ - t;,) (u,) - (v,) )(ix 

> ^\n, - v,f{u, + t;.)"$f„|(n, + v,)dx = C'^ j\ej,fp^,'^dx , 

with 

= 1 if a = -1 and = 1 + a if a G (-1,0) . 
The relative entropy inequality ()29p with /? = |a| becomes 



d I I 

— H\a,\[Ue,V^\f]+Caj pl'^jldx<l I j,{t,xfdx + Mi{v,ipm,-i) 

+\m\4f)\\L^H\^\[u,,v,\f] 

with 

M4(i/,(/p^,7) >2||at<&j„|(/)||i«=||<I>|,|(/)-<l>j„|(/)/||z.«. 

+ U^.^uml- + \\^Umf\\L- . 
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Integrating this inequality over the time interval [0, t] and applying Gron- 
wall's inequality and Proposition 13.11 as before, we obtain 

Ca f C Pl'^jldx 

Jo Jo 

Summarizing, we have proved 

Proposition 3.3. Let {u^,v^) be a solution of the scaled, generalized Carle- 
man system Q with admissible initial and boundary data. Then there exists 
a positive constant L = L{i^,ipm,a,t,u'^"',v^"') > such that 

r C p'^^^3l{t,x)dxdt<L, 
Jo Jo 

pT^3l{t,x)dxdt < L, 



10 Jo 

for each T > and e > if a ^ [Oi 1]; o,nd 

f-T fl 



[ [ pl^jUt,x)dxdt<L 
Jo Jo 



whenever a G [—1, 11. 



In the case a € [0, 1] , the second bound follows from the first and Propo- 
sition [3TT] by Holder's inequality. 

4. The nonlinear diffusion limit 

Let {u^"' , v^"") , (p~^ and (p~ be admissible initial and boundary data; then, 
for each e > 0, let (u^, v^) be the solution to the scaled, generalized Carleman 
system ([3|). 

It follows from Propositions 13. l l and 13.2] that, for each T > 0, one has 

IIPe|lL°°(o,r;L2(o,i)) ^ K {v, ifrn, Oi,T .u'"^ , v"') and 

l|je||L2([0,T]x[0,l]) ^ J{y,'fm,a:T,u''',v''^) 

for each e > 0. By the Banach-Alaoglu theorem, for each T > 

(31) the family /o^ is relatively compact in L°°(0, T; L^(0, 1)) weak-* 
while 

(32) the family is relatively compact in L^([0,T] x [0, 1]) weak. 
Summing both equations in Carleman's system implies that 

dtPe = -dxje 

so that 

dtpe is bounded in L'^{0,T; H-^{0, 1)) ; 
with the bound on coming from Proposition 13.21 this control implies that 

(33) pe is relatively compact in C{[0,T]; H^^ (0,1)) 
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by Arzela-Ascoli's theorem. 
Likewise, since 

Uf = and = , 

2 2 ' 

it follows that 

dx{Pe + eje) = -2p'^je " edt{pe + ej,) is bounded in L\[0, 1]; F"^(0, T)) 
d,{p, - ej,) = -2p'^j, + edtip, - eje) is bounded in L\[0, 1];H-\0, T)) 
because of Propositions I3.H 13.21 and 13. 31 Hence 

(34) p^ lb eje is relatively compact in C([0, 1]; H^^{0, T)) 

by Arzela-Ascoli's theorem. 

Because on the nonlinearities that are present both in the scaled Carle- 
man system ([3]) and in the limiting nonlinear diffusion equation ([6]), weak 
compactness results as above are not enough to pass to the limit as e — > 0. 
Strong compactness of the family p^ is obtained by the same argument 
as in [To] (see also [9]), which we recall below. 

Consider the vector fields 

Pe = (PeJe) and qe = {-pe, e^je) ■ 

By Propositions 13.11 and 13. 2^ both vector fields satisfy 

Pe and Qe are bounded in L^([0,T] x [0, 1];R^) • 

Summing the first two equations in the scaled Carleman system ([3]) shows 
that 

diVt,x Pe = dtPe + dxje = 

while 

curlt^a: qe = e^dtje + d^Pe = -^p'^je ■ 
By Proposition 13.31 we therefore have 

divj_2;Pe and curlt^^; bounded in L^([0,T] x [0, 1]) . 

Pick any sequence ^ such that 

- p and - j in L\% T] x [0, 1]) 

as n — > oo. By compensated compactness (the div-curl lemma in [14J), we 
find that 

2,2-2, 2 
Pen • 9en = -Per, + ^nJe„ P ' 1 = P 

in the sense of Radon measures on (0, T) x (0,1) as 0. Because of 

Proposition 13. H this implies that 

p1^ p^ in the sense of Radon measures 

which implies in turn that the family 

(35) Pe is relatively compact in L^([0, T] x [0, 1]) strong. 

Let then {p,j,q) be a weak limit point in L^([0,r] x [0, 1]) of the family 
{pe,je, pfje) as € — > — the cxisteuce of such limits points being guaranteed 
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by the bounds in Propositions 13.1^ 13.21 and 13.31 together with the Banach- 
Alaoglu Theorem. 

Because of (j33p and of the initial condition of the Carleman system ([3]), 
one has 

(36) peC{[0,T];H-\0,l)) and pj^^^ = n^" + w^" . 

Because of (j34p and of the boundary conditions of the scaled Carleman 
system ([3]), one has 

(37) p G C([0, l];H-H0, T)) and pj^^^ = 2^- , while = 2ip+ . 
Summing both equations in the scaled Carleman system shows that 

dtPe + dxje = 

so that, by passing to the limit in the relation above, one arrives at 

(38) dtp + 5,j = 

in the sense of distributions on R!j_ x (0, 1). 

Subtracting the second equation from the first in the scaled Carleman 
system shows that 

so that, by passing to the limit in the equation above, one finds that 

dxP = -2q . 

Assume first that a E [— liO]. By (|35p and the weak-strong continuity of 
the product of two functionsp 



so that 



J— 2P f^xP — 2{l-cx)^^P 



2F ^XH — 2{l-a) 

Substituting this in ([38|) . we see that p satisfies 

dtp-^dy~'' = 0, 

pit,l)=2ip+{t), 

p{0,x) =u''\x)+v'''{x). 

Since this problem has a unique solution, the whole sequence /o^ — > p in 
-L^([0,T] X [0, 1]) as e — > 0, by the compactness statement in ([35]) . 

The case a G (0, 1) is slightly more delicate. Since p > a.e. but 
may vanish, we set ?? > and, using as above the strong compactness 
statement ([35]) together with the weak-strong continuity of the product of 
two functions, 

1 _ • 

{p + rj)»'^~ {p + vY^ 



^ For / : R+ ^ R continuous and sublinear at infinity, if ^ p in and — ^ j in 
L\ f{Pe)U ^ /(P)J in -D'. 
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SO that 



We recall the argument in [9j. Letting i] ^ shows that 

On the other hand, the second bound in Proposition 13.31 implies that 



7' 

Jo 



^fdxdt < +00 



-j{t,z)dz 



so that 

j = a.e. wherever p = . 

Hence 

ip>oj = i 

so that the limiting nonlinear diffusion equation (j39p also holds in the case 
a e (0,1). 

It only remains to treat the case a = 1. Proceeding as above, we obtain 
instead the relation 

— 7— j = -^9^1n(r/ + p). 
Integrating this relation in x shows that 

nx 

ln(r? + p)(t,x) =ln(?? + v9+(t))-2 / „ 

so that, by Proposition 13.11 together with the bound 9?„^ < ip~^ < v, we see 
that 

I I sup I ln(r/ + x)|fitfix < +00 . 
JO Jo V 

Hence, by dominated convergence, Inp G -^^"'^([0, T] x [0, 1]) so that p > a.e. 
on [0,r] X [0,1]. Therefore 

^i^jin L2([0,r]x[0,l]) 

p + T] 

and 

\n{rj + p) Inp in L^([0,T] x [0,1]) 

as ?? — > 0, so that 

3 = -\dx\-a.p. 
Substituting this in ([H5|) . we see that p satisfies 

dtp — \d1 In p = , 
p{m=2^-{t), 

(40) 

p(t,l)=2v.+ (t), 

p(0,x) = n^"(x) + f^"(x) . 
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Again, since this problem has at most one solution, the whole sequence 
converges to that solution p in L^([0,T] x [0, 1]) because of the compactness 
statement (|35]) . 

Thus we have established the nonlinear diffusion limit for each a € [— 1, 1]. 
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